Beam-Breakup Instability Theory for Energy Recovery Linacs 
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Here we will derive the general theory of the beam-breakup instability in recirculating linear 
accelerators, in which the bunches do not have to be at the same RF phase during each recirculation 
turn. This is important for the description of energy recovery linacs (ERLs) where bunches are 
recirculated at a decelerating phase of the RF wave and for other recirculator arrangements where 
different RF phases are of an advantage. Furthermore it can be used for the analysis of phase errors 
of recirculated bunches. It is shown how the threshold current for a given linac can be computed 
and a remarkable agreement with tracking data is demonstrated. The general formulas are then 
analyzed for several analytically solvable cases, which show: (a) Why different higher order modes 
(HOM) in one cavity do not couple so that the most dangerous modes can be considered individually, 
(b) How different HOM frequencies have to be in order to consider them separately, (c) That no 
optics can cause the HOMs of two cavities to cancel, (d) How an optics can avoid the addition of 
the instabilities of two cavities, (e) How a HOM in a multiple-turn recirculator interferes with itself. 
Furthermore, a simple method to compute the orbit deviations produced by cavity misalignments 
has also been introduced. It is shown that the BBU instability always occurs before the orbit 
excursion becomes very large. 



I. INTRODUCTION 

Synchrotron light sources based on Energy Recovery 
Linacs (ERLs) show promise to deliver X-ray beams with 
both brilliance and X-ray pulse duration far superior to 
the values that can be achieved with storage ring technol- 
ogy. This is due to the fact that the emittances in an ERL 
are largely determined by a laser-driven source, technol- 
ogy which has been improving steadily over the years, 
and will undoubtedly improve further. To generate high 
brilliance high flux X-rays it is necessary to accelerate 
beams to the energies (several GeV) and with the cur- 
rents (several 100 mA) that are typical in these storage 
rings. This would require that the linac delivers a power 
of order 1 GW to the beam. Without somehow recovering 
this energy after the beam has been used, such a device 
would be practically unfeasible. Energy recovery yij can 
be achieved by decelerating high energy electrons to gen- 
erate cavity fields which in turn accelerate new electrons 
to high energy. With this, large beam powers that are 
not accessible in a conventional linac, can be produced. 

Several laboratories have proposed high power ERLs 
for different purposes. Designs for light production with 
different parameter sets and various applications are be- 
ing worked on by Cornell University 0, > BNL 
Daresbury [I, TJNAF @, JAERI 0, the University of 
Erlangen [J; Novosibirsk @, and KEK [nj. TJNAF 
has incorporated an ERL in its design of an clcctron- 
ion collider (EIC) [ll| for medium energy physics, while 
BNL is working on an ERL-based electron cooler [lj for 
the ions in the relativistic ion collider (RHIC). The work 
at TJNAF, JAERI and Novosibirsk is based on existing 
ERLs of relatively small scale. 

One important limitation to the current that can be ac- 
celerated in such an ERL is given by the beam-breakup 
(BBU) instability. The size and cost of all these new ac- 
celerators certainly requires a very detailed understand- 
ing of this limitation. 



A theory of BBU instability in recirculating linacs, 
where the energy is not recovered in the linac, but where 
energy is added to the beam when it returns after each 
recirculation turn, was presented in [13|. Such a linac 
can consists of many cavities and several recirculation 
turns can be used. This original theory was additionally 
restricted to scenarios where the bunches of the different 
turns are in the linac at about the same accelerating RF 
phase, such as in the so-called continuous wave (CW) op- 
eration where every bucket is filled. Tracking simulations 
[l4| compared well with this theory. In the following we 
therefore refer to it as the CW recirculator BBU theory. 
It determines above what threshold current ith the trans- 
verse bunch position x displays undamped oscillations in 
the presence of a higher order mode (HOM) with fre- 
quency u\ . If there is only one higher order mode and one 
recirculation turn with a recirculation time t r in the linac, 
the following formula is obtained for Tyi sinuj\t r < 0: 



ith = - 



2c 2 
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(1) 



where c is the speed of light, T 12 is the element of the 
transport matrix that relates initial transverse momen- 
tum p x before and x after the recirculation loop, e is the 
elementary charge, (R/Q)\Q\ is the impedance (in units 
of f2) of the higher order mode driving the instability, 
Q\ is its quality factor. A corresponding formula had 
already been presented in [2]]]. Occasionally, additional 
factors are found when this equation is stated [Tfil flit IT^ | . 
We give a concise derivation which shows that no such 
additional factors are required. 

In the past efforts have been made to derive threshold 
currents from the analysis of experimental data obtained 
at the TJNAF FEL-ERL OEai using the CW recircula- 
tor BBU theory. The data has not been interpreted satis- 
factorily, and at least part of the reason could be that the 
underlying theory had not been derived for ERL opera- 
tion. The theory presented in this paper should therefore 
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be used to extend and improve these previous analysis re- 
sults, as it is directly applicable to ERL operation. 

Here we will derive the general theory of the beam- 
breakup instability in recirculating linear accelerators, in 
which the bunches do not have to be at the same RF 
phase during each recirculation turn, similar to what has 
first been presented in [20|. First we treat the simplest 
case of one dipole HOM and one recirculation turn. Then 
we allow many HOMs and many recirculations, and fi- 
nally we analyze several analytically solvable cases. 

II. ONE DIPOLE HOM AND ONE 
RECIRCULATION 

For recirculating linacs, the simplest case of one HOM 
and one recirculation loop has long been described [l7| . 
Previous theories which assume that the recirculation 
time is an integral number of RF periods should not be 
used when investigating ERLs. Next we derive more gen- 
eral formulas that may be applied for arbitrary recircu- 
lation times. 



A. The Dispersion Relation 

In the simplest model of multi-pass beam breakup, 
bunches are injected into a cavity, which is assumed to 
have one dipole HOM (e.g. TMn-like mode), accelerated 
in the cavity and then recirculated to pass the cavity a 
second time before they are ejected. In the case of a two- 
turn recirculating linac, each bunch would be accelerated 
on both passes through the linac and ejected to a user 
area. The RF phase of the bunch would therefore be ap- 
proximately the same on both passes. In an ERL, the RF 
phase on the second pass through the cavity is shifted by 
7r with respect to the first pass so that the energy that 
the bunch gains in the first pass is returned to the cav- 
ity during the second pass and the bunch is ejected with 
reduced energy into a beam dump. 

If a dipole HOM is excited in the cavity, then a bunch 
that enters the cavity on axis experiences a transverse 
kick and starts to oscillate around the design orbit of the 
recirculation loop and returns to the cavity with a trans- 
verse offset. This offset leads to a change in the energy 
of the HOM. If it increases the HOM energy, transverse 
kicks experienced by subsequent bunches will be larger, 
which will in turn lead to a further growth of the HOM 
energy once the kicked bunch returns to the linac: an 
instability develops. 

To describe this effect, we use the ideas and nomen- 
clature from 0. When a current I(t') passes the cavity 
during its recovery loop at a time t', the charge I(t')dt' 
with the transverse offset x(t') excites the dipole HOM, 
creating a transverse momentum for particles traveling 
through the cavity subsequently at time t, 

Ap x (t) = -W(t - t')x{t')I{t')dt' , (2) 



where the wake function W(t) describes the transverse 
force at time r after the HOM was excited. The mo- 
mentum transfer is described by an effective transverse 
voltage of the HOM, V{t) = § Ap x (t). 

Assuming that all bunches are injected on the cav- 
ity's central axis, they do not excite dipole HOMs on 
their first path through the cavity. However, the effec- 
tive transverse voltage of the HOM determines what kick 
Ap x (t) the bunch sees and what position it will have 
when it returns to the cavity after the recirculation time 
t r . The transfer matrix element T 12 maps the transverse 
momentum p x (t) to x(t + t r ) = T\2p x {t). Inserting this 
into Eq. (J2J) leads to an integral equation for the HOM's 
effective voltage, 

V(t) = / W(t - t')I(t')T 12 -V(t' - t r )dt' . (3) 

J-oo c 

To solve this integral equation, one assumes that the 
current is a continuous stream of short pulses being in- 
jected at multiples of an interval between bunches tb, so 
that the current on the second turn is given by 

oo 

I(t)=I t b S D (t - t r - mt b ) , (4) 

ra— — oo 

(5d being the Dirac-delta function. Note that i& is an 
integer multiple of the RF circulation time to = 2tt/luo 
for the RF frequency luq. We write the recirculation time 
in terms of the time tb between bunches as 

U = (n r - S)t b , (5) 

with an integer n r and 5 £ [0,1). For a recirculating linac 
one has 

Stb « nto , (6) 

and for an ERL one has Stb ~ (n + \)t$ for some integer 
n. A "+" sign in Eq. 0that defines 8 may seem more 
natural but our choice leads to simplified equations. 

The HOM voltage at a time t S [ntb + t r ,ntb + t r + tb) 
is given by 

n 

V(t) = I t b T 12 - Y W{t-t r -mtb)V(mt b ) . (7) 
c * — ' 

m= — oo 

Evaluating this at the time t = ntb + t r when the recir- 
culated bunches pass through the cavity leads to 

oo 

Vintb + t r ) = I t b T 12 - V W(mt b )V([n - m]t b ) . (8) 
c z — ' 

m— 

In the CW recirculator BBU theory this difference equa- 
tion was dealt with by assuming that the voltage can be 
written as V(t) — V~oe~ luJt for t — ntb where a positive 
imaginary part of u> indicates instability. Note that this 
does not require V(t) to be a harmonic function, but that 
it can be a linear combination of harmonic functions with 
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frequencies us + ra^- for integers m. This distinction has 
not been always made clear and is a potential source of 
confusion jl^ ■ One obtains the equation 



~T = t b li2-e > 

in C '-^ 

u m=0 



W{mt b )e tujmtb 



(9) 



The smallest value of the current Iq for which there is a 
real u> is the threshold current 7th of the instability. 

We proceed by writing V(t) in terms of its Laplace 
transform, retaining all possible frequencies in HOM volt- 
age, which automatically enables proper description of 
arbitrary recirculating configuration: 



V{t) = h 

ITT 



V(uj')e- luJ,t duj' . 



(10) 



Note that this is not the conventional way of writing a 
Laplace transform. We have chosen this notation in order 
to make it appear more similar to a Fourier Transform. 
It also makes the subsequent notation more similar to the 
CW recirculator BBU theory. The Laplace transform is 
used rather than the Fourier transform since we want to 
analyze the onset of instability where the frequencies u 
become complex. With the following definition 



E 



2tt 
— 1 

t b 



we obtain 



V(nt b )e 



iujntb 



(ii) 



(12) 



Since V s (w) is periodic with 2ir/t b , it has a Fourier series, 
and its Fourier coefficients are V(ntb). This shows that 
V s (ui) does not vanish. We can therefore choose t = 
(n + n r )tb in Eq. J7J) and sum over n, 



n+n r ]tb 



(13) 



A)£fcTi 2 - 

c 

oo oo 



x J2 W([m + S}t b )V([n~m}t b )e l ^ n+n ^ tb 

n— — oo m— 

oo 

= I t b T 12 -e lujnrtb V W([m + S}t b )e^ mt "V^(uj) . 

r — ' 



This finally yields the dispersion relation between Iq and 
u> which can be used for all 8. A corresponding derivation, 
which has treated the beam recirculation in a way that 
can be applied to ERLs, has been presented in [2(|. We 
believe that this paper should be referenced more often, 
since it is hardly referenced, whereas the earlier papers 
with CW recirculator BBU theory, which was not derived 
for ERLs, are often referenced in the context of ERLs, 
where they are not strictly applicable. 



The Far-Field Wake 



The sum in the dispersion relation 



i = t b T 12 -e Mb w(8) , 
I c 



j(5) = J2W{[n + S]t b )e 



(14) 



(15) 



can be computed when the far-field approximation for 
the wake function is used, 



W ( T ) = ( 7;) ^ A e"^t T sinw A r . (16) 



QJx 2c' 

With ujf — u\ ± i^Q^ an d lo + ~ lo + ^he recm i re( i 
sum can be evaluated if lm(uj) > and becomes 

oo 

(6) = ^W{[n + 8\t b )e™ nt » (17) 



w 



4ic 



~St h 



QJx 4c 

e iw+(5-i)t b S j n (u x 8t b ) - e lSuJ+t » sin(w A [£ - l\t b ) 
x — 

COSW+ifc — COSLO\t b 

The dispersion relation thus becomes 



a Stb [cos(uj + t b ) — cos(uj\t b )] 



Jo 



-iu)n r tb 



JCT 12 



e -^ +t » sm(Su\t b ) - sin([<5 - l]iu x t b ) ' 

(18) 

with K, = *b^(0)A^2 L - For 8 — this becomes the dis- 
persion relation of the CW recirculator BBU theory: 



e -iut r cos(uj + t b ) - cos(L0 X t b ) , 1Q . 



K,T X 



sin(w A t h ) 



This describes the case when the recirculated bunches are 
in the same buckets as the accelerated bunches. When 
the recirculated bunches are just between accelerated 
bunches, then 8 = ^, 



Jo 



12 



_ e -iut r cos(uj + t b ) - cos(uxt b ) 



cos(oj + t b /2) sm(uj\t b /2) 



For the case that every bucket is filled, this would be an 
ERL with t b = to. The dispersion relation for ERLs for 
other 8 is less simple than Eq. 120fl and has 8t b — (n + 
i)^o in Eq. I|18fl . This occurs when the decelerating and 
accelerating bunches are not perfectly centered between 
each other. 
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FIG. 2: Awt{, (blue dots) and ecotu;t r (blue and green 
curves). The curve is dark blue in the region where 
T\2 sin Lot r < and light green otherwise. The to of the in- 
stability is given by the intersection of the dotted line with a 
dark blue solid curve for which Aoj is smallest. Left: A case 
of T12 smto\t r < 0. Right: A case of T12 smco\t r > 0. 



FIG. 1: Iq(io) in the complex plain for u £ [0,ir/tb]. The scale 
is arbitrary. 



C. The Threshold Current 

For a given positive current Iq, the values of lo that 
satisfy the dispersion relation Eq. I |18(l will in general 
be complex. If they all have negative imaginary parts, 
the beam motion is stable. If one of them has positive 
imaginary part it will be unstable. 

For small currents the beam motion is stable. When 
the current is increased, at some point, one of these co 
will become real. At this point the threshold current 
is reached. The threshold current ith is therefore the 
smallest current Iq for which there is a real lo that satisfies 
the dispersion relation. To find this current, we note that 

I (w + ^) = JoH , Io(-w*) = I » . (21) 

It is therefore sufficient to investigate co <E [0,7r/tf,]. 

Figure ^ shows Iq(lo) in the complex plain for co £ 
[0, Tr/tb]. The intersection with the real axis that has the 
smallest positive value yields the threshold current. 



III. APPROXIMATE THRESHOLD CURRENT 

It can often be justified to linearize in e = ^j^tb, e <C 1, 
which describes a situation when HOM decay is negli- 
gible on the time scale of the bunch spacing tb- This 
applies to linacs when nearly every RF bucket is filled. 
The smallest |/o| m Eq. H18fl is obtained when cosojtb 
is close to cosui\tb, which occurs whenever w is close to 



Y~ for any integer n. Due to Eq. (|2*T|l . 

all these frequencies lead to the same threshold current. 
We therefore additionally linearize in Aw = lo — uj\. n± , 



assuming Aojtb -€i 1. This leads to 
2 



Io = T 



JCT 12 



- lun - tb e^ x - 0±Stb {Aujtb + ie) 



(22) 



Within the linearization, the phase factors could be com- 



bined to e 1 ' 



This is not done in order to retain 



the symmetries of Eq. (|21[) for uj — > to + 2ir/tb and for 
lu —> —lo*, i.e. Alo — ► —Alo* and L0x, n ± ~> ^>X,—n^- 

Due to these symmetries, the real current close to 
LO\^ n ± is the same for each of these frequencies. With- 
out loss of generality we therefore use to\. n ± — lo\ and 
no longer require the symmetries, 



Io = - 



12 



'(Acotb + it) 



(23) 



Since Iq must be real, Acotb sin iot r = e cos cot r leading to 
the following two equivalent equations, 



K T12 sin tot, 
lo = —sign (sin u>t r ) 



(AloU) 2 



(24) 
(25) 



For this formula to describe the threshold current, it is 
required that lo > and therefore T12 smcot r < 0. 

Figure [3 shows Acotb and ecotcot r versus cotb G [0,tt] 
for t r — 6.88ifc. The dotted line and the curve have to 
meet at a region where T±% sinw^t < 0, which is indicated 
by a dark blue curve. 

When n r e <C 1, i.e. HOM decay is negligible on the 
time scale of recirculation time, then co ss co\ in the region 
where T12 sm.co\t r < and one obtains 



/n = 



K Tn sa\co\t r 



(26) 



which is the traditional and commonly used approxima- 
tion in Eq. Q which had been derived for 6 = 0. 



5 



ln(I th [A]) 




FIG. 3: Threshold current obtained by tracking (red dots) 
and approximate analytical solution (top) and by a numerical 
solution (bottom) of the dispersion relation Eq. 1181 . Param- 
eters: rir-S € [6.135,7.234], (R/Q)\ = 100 SI, Q\ = 10 4 , 
T12 = -l(T 6 eV/c, coxtt = 9.67. 



In the region with Tyi %va.oo\t r > one obtains ton r 
nir, which can be used in Eq. I|25|) . 




— Mod(oj\t r , 7r) 
nt 



(27) 



For n r e ^> 1, i.e. when HOM damping is substantial on 
the recirculation time scale, one again uses Eq. 1251) with 
uit r ss (2n=F i)-7r, where the + or — sign is determined by 
the sign of T 12 smu>\t r , 




(28) 



Note that in this case the threshold current weakly de- 
pends on t r and can be estimated simply by Iq = 
2e/lC\T 12 \. 

One can perform these approximations more accu- 
rately, for example, by approximating cot cot r by a line 
or second order curve. However, in regions where n r e is 
not much larger or much smaller than 1, simple formulas 
cannot be found. 

Figure |3| (top) shows the threshold current obtained 
with the approximate analytic solution compared with 
the threshold current that is found by tracking particles 
for the simple case of one cavity with one HOM and one 
recirculation loop. 



Figure (bottom) compares the same tracking re- 
sults with a numerical solution of the dispersion rela- 
tion Eq. I|18l) . The data agrees remarkably well with the 
approximate formula in the region where Alo is small, 
i.e. where Jo is relatively small. In the region where the 
threshold current is relatively large, the agreement with 
the approximate expression is not satisfactory, however. 

To find the threshold current with Eq. I|18|) , the small- 
est positive real value of Jo for to e [0, 7r/tf,] was found 
by linearly interpolating 100 points in the region to G 

K - — r - h ,u x + —-\. 



A. Instability Growth Rate 

We denote the threshold current by Jo and the real 
frequency Co 6 [0, 7r/t b ] satisfies Eq. ljT%|) for this current. 
When the current Jo is slightly larger than Jo, there is 
one frequency a; (Jo) that satisfies Eq. (|18f) and is close to 
Co. It has a positive imaginary part. All other frequen- 
cies lo at which Eq. (|18fl holds and for which therefore 
V s (o>) might not vanish have an imaginary part that is 
not positive. In Eq. (|10|) there are therefore exponen- 
tially growing terms. For currents that are only slightly 
larger than the threshold current, the complex frequency 
lo(Iq) can be expanded with respect to A J = Jq — Jo, 



(29) 



The rise time per current of the instability is thus given 
by a — 3{^}|/ Q • The dispersion relation Eq. l(T%|) leads 
to a long formula. However, using the simplified Eq. i|22|) 
leads to 



a 



1 



1 



d Jo / dco 



(30) 



T 3{( 
Jo 

1 



in r tt. 



4c 



{Cb - LO\)tb 

(J /CTi 2 )V 



Jo^b 4 + 8?i r e + (Jo/CTi 2 n r ) 2 

Provided parameters are not in the region where the 
curves diverge in Fig. El (top), i.e. sinw>< r is not close 
to zero, the following approximate formulas hold for the 
growth rate: For n r e ^ 1 one obtains a — -^—e and for 

lotb 

n r e ^> 1 one obtains a = — —. 



IV. MULTIPLE DIPOLE HOMS AND 
MULTIPLE RECIRCULATIONS 

Recirculating linacs with many cavities and several re- 
circulation loops have been considered early on 0, 0] . 
Here we use the same nomenclature as much as possible. 
The N higher order modes, which can be associated with 
different cavities, are numbered by an index i. The N p 
passes through the linac are numbered by an index J. 
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The horizontal position and momentum that the beam 
has at time t in the HOM i during turn I is denoted 
z\(t) — (x\ (t),pj (t)). The transport matrix that trans- 
ports the phase space vector z 1 - at HOM j during turn J 

to z\ is denoted and the time it takes to transport 
a particle from the beginning of the first turn to HOM 
i during turn / is denoted t\. The beam is propagated 
from after HOM % - 1 to after HOM % by 



^(*) = T^_ 1 -2f_i(*-[*i - 



(31) 



This equation can be iterated to obtain the phase space 
coordinates as a function of the HOM strength that cre- 
ates the orbit oscillations. With the matrix element 
T-J = (Tf/)i 2 one obtains 



7-1 N 

7=1 3 =1 ° 



The strength Vi (t) of the HOM i is created by all parti- 
cles that have traveled through that HOM via the integral 



Vi(t) = £Wi(t- Otf . (33) 

J-oo I=1 

where // (t) is the current at time t that the fraction of 
the beam has which passes the HOM i on turn /. Com- 
bining this with Eq. (|32() leads to the following integral- 
difference equation: 



J — OO j i 



(34) 



7=1 

I N L ,{i-l) 



N u (i - 1) 



' 7=1 3=1 

fiV, if/^J; 



(35) 



Now the approximation of short bunches is used. The 
current is given at time i by pulses that are equally spaced 
with the distance t b , 



1} (t) = E IotbS ^ mtfc ) 



(36) 



m= — oo 



This reduces the integral to a sum, 

n{t,t{) N p 

Vi(t) = -I t b J^Wiit-tl-mh) (37) 



m— — oo 1=1 
I Njj(i-1) 
x E E TffVj(mt b + t 

J=l j=l 



where n(t, tj) — Max m (t > mt b + tj). Computing 

oo 

V t L '= E ^i(n*fc + *f)e Wb (38) 

n= — oo 

leads to 



oo oo 



^ = E E E w ^ mt * + ^ - ** ) (3 9 ) 

n=— oo m=m* 7=1 
/ iVjj(i-l) 

X E E ^ ([n-m]t b + i/)e^ . 

J=l 3=1 

The second summation starts at m* = n — n(nt b + , ) 
which can be simplified by writing t\ = (nf — Sf)t b , 

m* = n-n([n + nf -nj]f 6 ,[5r^]4) (40) 
= -Max TO ([nf - nf]t 6 > [m + 5f - 5{]t b ) 
= nl-nf + ii(I,L) , 

where ji(I,L) = 1 if 8f > Sj and otherwise. Shifting 
the summation index m now leads to 



co oo Np 



v " = -J^ E E E (4D 

n= — oo m=0 1=1 

I N L ,(i-l) 

^am+ 7i (/,i)+«5/-^6)E E T " 

J=l 3 = 1 

x V;j{\n-m- n( + nf- 7i (J, L)]i 6 + tj )e lujntb . 

and with 5i(I, L) = L) + Sf — 5f , which is between 
and 1, shifting the index n finally leads to the relation 



oo Np 



vt = -ttEE^d™ 1 ^ 1 '*) ( 42 ) 

m=0 7=1 

x jT/J e iai[m+nJ-ref+7 4 (/,L)]t»yJ 
J=l 3=1 

The following sum is equivalent to that in Eq. i|17|) . 



,(*) = E Wi([m + %)( 



(43) 



m=0 



Equation (|42|) reduces to 



1^ = -t 6 Vm J (<5 l (/,L)) e i 

/0 C 7=1 

7 Atjj(i-l) 

x E E W ■ 

7=1 3=1 



(44) 



If a vector V is introduced that has the coefficients V/ , 
this equation can be written in matrix form, 



-y = w( w )uy , 



(45) 
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where the matrix M = W(w)U is determined by 
Eq. H44f> . When all electrons are considered to have the 
speed of light, t\ — tf does not depend on the HOM 
number i and we therefore drop this index and obtain 
the following matrix coefficients: 

n p 

Mi y = U b J2 Wi(S(I,L))e iuTop ^ )U T^ 



0, otherwise, 



(46) 



where Top(a;) is the smallest integer that is equal to or 
larger than x. With Kronecker 6ik this determines the 
matrices W and U to be 



Wk 1 



•i hWl (<5(/,L))e MTop( ^ )t ^ lfc , (47) 



U, 



u 



(48) 



For each frequency to, Iq is an eigenvalue of M(w). 
Since the eigenvalues are in general complex, but Iq has 
to be real, the threshold current is determined by the 
largest real eigenvalue of M.(uo). The matrix has the 
properties 



M(co- 



2tt, 



= M(w) , M(-w*) = M*(w) 



(49) 



and it is therefore again sufficient to investigate lo 6 
[0,7r/tf,] to find the threshold current. 



Note that V^ F never appears on the right hand side of 
Eq. I)42|l so that the dimension of M can be reduced by 
one to N x N p — 1. Furthermore the dimension can be 
reduced when two fractional parts 6( and 6j are equal 
since then V/ and Vj 7 are identical. Note also that for 
N = 1 and N p — 2 Eq. I|42[l reduces to the dispersion 
relation for one HOM in Eq. (JTHJ) . 



A. Instability Growth Rate 

The growth rate of the instability is again computed 
by first obtaining the threshold current io and the real 
frequency lo for which 1$ 1 is an eigenvalue of M(w). If 
this is the fcth eigenvalue Afc(w) of the matrix M(u), then 
the growth rate of the instability is given by 



a = —Ir 



-2 cv 



duo 



(50) 



V. MULTIPLE HOMS IN ONE CAVITY 

The presented theory for multiple HOMs and multi- 
ple recirculation turns can in general only be evaluated 
with computers. However, for some simple situation an 
analytical understanding is possible. 



One such case is an accelerator with one recirculation 
loop and one cavity in which many HOMs can be ex- 
cited. The (1,2) matrix elements that refer to transport 
between HOMs for the same pass are zero, T^ 1 = 0. All 
matrix elements that describe the recirculation loop are 
identical, T?j — Tvz i£ I ^ J. The matrix elements in 
Eq. I|46fl are then given by 



M t f = ,M t f = - c t bWi (S(2, i)) e ia,Top(£ ^* 6) T 12 



(51) 

Equation (|45|) becomes with 6 — 6(2, 1), 

T V ^ = wifie^^TnV} • (52) 

io c t-r 1 J 

3 

A summation over the index i shows that • can only 
be nonzero when 



— = thTvx-e 
io c 



iajn r tt, 



Wi(6) . 



(53) 



Comparing this with Eq. (|14fl shows that one only has to 
replace w(6) of the single HOM by J2i w i($) to arrive at 
the threshold formula for the multi-HOM case. To find 
the smallest real 7o that Eq. I|53|l can produce for real 
lo, Wi(5) has to be maximized. The sum is especially 
large when the denominator of one of the terms is very 
small, i.e. when coswif, « cos LO\ti,. 

When the HOM frequencies modulo 2n/tb are suffi- 
ciently different for the different HOMs, the maximal ab- 
solute value of J2i w i(d) wm be close to the largest abso- 
lute value that any of the Wi (6) could have individually. 
This is due to the fact that all the Wj(6) for j ^ i are rel- 
atively small for frequencies to for which the denominator 
of Wi(5) is small. 

For HOM frequencies that are sufficiently different in 
the above sense, the threshold current for several HOMs 
therefore does not differ significantly from the threshold 
current of the worst individual HOM. 

Figure shows how the threshold current changes 
when the frequency of one HOM is fixed at a small thresh- 
old current with |sina>it r | = 1 and a second HOM fre- 
quency is varied. Superimposed is I t h if only the sec- 
ond HOM is present. It is apparent that the HOM that 
would produce the larger threshold-current if it was solely 
present only influences the threshold current of the pair 
when the frequencies ±wi mod27r/tb and ±002 mod27r /% 
are closer together than about Aco\ = e/tb = lo\/2Q\. 

One can draw the conclusion that in the case of 
many HOMs they do not interact destructively when the 
frequencies ±00 \ mod 2ujt\, are not very close together. 
Tracking simulations also demonstrate this effect. 

One strategy to increase the BBU threshold current 
is the introduction of HOM frequency spreads between 
cavities. As an example, Fig. [5] shows the threshold cur- 
rent found by tracking as a function of uniform frequency 
spread of 20 HOMs in a single cavity. For all three curves 
(R/Q)\Q\ is the same. It is seen that for lower Q\, the 
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Wth[A]) 





FIG. 4: Dark black curve: the threshold current ln(/ t h[A]) for 
one HOM at u>i /2n — 2 GHz as a function of a second HOM 
with frequency ui2- Light green curve: threshold current when 
only the second HOM is present. Light red lines: frequencies 
for which cosu^ij ~ cosoiit;, where the threshold current is 
not simply the minimum of the threshold currents produced 
by the individual HOMs. 




frequency spread (MHz) 

FIG. 5: The threshold current Ith as a function of uniform 
equidistant frequency spread for 20 HOMs in a single cav- 
ity. Abscissa displays frequency difference between the two 
HOMs adjacent in frequency. Parameters: lu\/2tt = 2 GHz, 
(R/Q)xQx = 5000 ST2, T 12 = -10 _6 eV/c, t r = 1.000125 • 
10 -6 s, uiq/2-k = 1.3 GHz. Threshold is determined by track- 
ing with accuracy 0.1%. 



curve begins to saturate for larger frequency spread than 
for the high Q\ case. It is also seen that the threshold 
for frequency spread uj\/2Q\ is similar in all three cases. 
Both observations are consistent with the above assertion 
that HOMs do not interfere when they are further apart 
than oj\/2Q\. The fact that oscillations in the Fig. [5] 
are smaller for low Q\ is consistent with the conclusion 
that a broader HOM resonance peak should lead to more 



overlaps and as a result to less pronounced differences in 
the threshold for different HOM frequencies. 



VI. ONE DIPOLE HOM IN TWO CAVITIES 

In order to see how two dipole HOMs interact, we will 
now analyze a set of two HOMs with one recirculation 
loop, i.e. N p — 2. We abbreviate <5 = 5(2,1) and a r = 
e zum r t b _ rp^g di S p ers i on relation is then given by 



= -4N 

c 



V, 1 



(54) 



u>2 [GH Z ] where the matrix N is given by 



Wl (5)a r T$ 
w 2 (0)Ti} + w 2 (5)a r T$l 



wi{5)a r T$ 
W2(S)a r T22 



w 2 (5)a r Ti? 




(55) 

The third row is similar to the first row and eliminating 
it by similarity transformations leads to the 2x2 matrix 



w\{S)a r 







"2(0) 



io 2 (<5) 



U T 11 



w 2 (5)a r 

rp2\ 
± \\ 
rp 2 \ 
J 21 



(56) 



Wl(0) y22 
wi(5)a r 21 



rp21 

1 12 

rp21 

1 22 



Even though u>i(5) appears in the denominator, the 
formula for the eigenvalues of this matrix does not con- 
tain such a denominator. Therefore the largest real eigen- 
value will again occur at a frequency u> for which one 
of the HOM frequencies satisfies cos(u;tb) f» cos(uj\tb). 
The term Wi(6) and tUj(0) of the other HOM can then 
again be neglected, so that for sufficiently different HOM 
frequencies mod 27r/i the threshold current is again ap- 
proximately determined by the HOM which would have 
the smallest Ith if there were no other HOMs present. 

We therefore now assume that the two HOMs are 
equal, wi = u>2- Now one can perform an approxima- 
tion analogous to Eq. (|2*5|) leading to 



N 



K, 



1 



rp21 

1 21 



21 



rp2\ 

\ 1 21 



(57) 



rp21 

1 12 

7^22-1 rp21 
1 21> ± 22 



It is interesting to analyze whether the effect of the 
two HOMs can cancel. A cancellation could occur most 
naturally when the linac and the recovery loop are mirror 
symmetric. A mirror symmetry of the linac means that 
the beta functions of the first pass, going from low to 
high energy, are the mirror image of those of the second 
pass, going from high to low energy. An example of such 
an optics is shown in Fig. 

Since the (1,2) element of the transport matrix be- 
tween a region with momentum po and a region with 
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momentum p can be written with Twiss parameters as 



T U = J ^° sin A* 
V PPo 



(58) 



this symmetry leads to T 22 = T^j . An additional mirror 
symmetry of the return arc leads to = T\\- The 
eigenvalues of the matrix in Eq. I|57(l then become, 



1 



.—pi"' 

2 

T 21 4- 



(59) 



2e 



-iwt r T^l^ 
J 21 J 



Since there are two solutions to the quadratic eigenvalue 
equation, to every eigenvalue that is smaller than 1/ith 
of a single cavity, there exists in general the one that is 
larger. Therefore two cavities do not compensate their 
instabilities, but it is possible to decouple the cavities to 
the extent that the combined threshold current is just 
as large as that for a single cavity. For this, one has to 
choose rf 2 = 0, i.e. the phase advance of the return arc 
has to be a multiple of tt. 

Since the kick of a HOM disturbs the beam most at 
low energy, the first and the last cavity of an ERL are the 
strongest contributors to BBU. It seems therefore advis- 
able to adjust the phase advance of the arc to a multiple 
of 7r also when the linac has more than two cavities. 



A. Multiple Recirculation Turns 

One could envision a multi-turn recirculating linac as 
ERL. The beam would pass the same linac N r times to 
reach its top energy, and subsequently it would be decel- 
erated in just as many turns through the linac. The cur- 
rent in the cavities would be 2N r times higher than the 
current that is available at high energy. One could there- 
fore conjecture that the BBU threshold current is N r 
times smaller than for a one-turn ERL. Here we will show 
that the threshold current can be significantly smaller 
than that conjecture and in general can be expected to 
decrease quadratically with N r . 




150 200 250 
position in linac (m) 



FIG. 6: Example of a mirror-symmetric linac optics. Dark red 
curve (light green curve) is horizontal (vertical) beta function 
in the linac for accelerating beam. Dashed lines show lattice 
function for decelerating beam. 



Each bunch passes the linac N p = 2N r times and the 
matrix in Eq. (|46[l becomes 



M 



LJ 



-t b w(S(I,L))e 2ujTop{ ^ )tb T IJ , (60) 



7=7+1 



where the lower indexes have been suppressed since only 
one HOM is considered. The eigenvalue equation thus 
becomes 



N B N„ 



±-V L = - c t b J2 ]T W ( ( 5(/,L))e^ Top( ^ ) * 6 T /J F J 



j=i i=j+i 



(61) 

An approximation equivalent to that leading to 
Eqs. (HHJl and (^3 results in 



Io 2 Aujtt, + ie 



(62) 



N p N p 

E E - 

,7=1 7=J+1 



J [Top(i^^)-6(I,L)]tb T UyJ 



Comparing the exponent on the right hand side to those 
in Eqs. (|44l) and 1461) shows that this can be written as 



K 



1 



— V L e iwtL — 

Io 2 Autb + it 

Np Np 

e iu(t I -t J ) T u v 

J=l 7=7+1 



(63) 



Since the right hand side does not depend on L, all the 
terms V J e luJt are equivalent and the condition that they 
do not vanish is 



/C 1 



1 



Jo 2 Aujti, + ie 



Ju(t' -t J )rpIJ 



T , (64) 



7>7 



where the double sum goes over all pairs of / and J for 
which I > J. Similar to the condition obtained from 
Eq. (|23(l which is analyzed with Fig. the fact that Iq 
has to be real entails the condition 



7>7 



cosH^ ~t J ])T IJ 



7>7 



_ Alu 

sin(a;[t' ~t J ])T" ~ ~ 



(65) 



In regions where J^7>7 sin(o-'[i / — t J ])T IJ > 0, we again 
obtain the approximation that to ~ lo\. The equation for 
the threshold current of an N r times recirculating ERL 



corresponds therefore to that of the case without 
recirculation in Eq. 

2c 2 1 



jN r _ _ 
x th ~ 



e(f)AQA^A E sinM^-F])T" 



(66) 



7>7 



A comparison with Eq. Q shows that this current 
is smaller than the one-turn ERL by a factor of up to 
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5^ 7>J iT^I/lTjd- This is in agreement with earlier result 
presented in |2J . Assuming that all matrix elements are 
of about equal magnitude, the threshold current in an N r 
times recirculating ERL is therefore in general smaller by 
about a factor of N r (2N r — 1). This conclusion is con- 
sistent with tracking results for microtrons (2^ an d for 
two-turn ERL [2*3 . The scaling in a particular case, how- 
ever, can be quite different depending on details of the 
lattice design, e.g. approximate scaling with N r was re- 
ported in |24| . 



VII. CAVITY MISALIGNMENTS 

In the derivation above, it was assumed that the 
bunches travel along the cavities' symmetry axes when 
the current is below the threshold for BBU instability. 
When the cavities are misaligned, the beam will excite 
dipole higher order modes even below the threshold and 
the trajectory will be disturbed by these modes. 

Let us assume that the ith cavity is misaligned with 
respect to the path adjustment of the Jth turn by x$ i: 
leading to the misalignment vector xq . The dipole HOMs 
that are excited by the beam are now not only due to the 
beam position fluctuation that is produced by the HOMs 
themselves, but additionally due to the cavity misalign- 
ments. The HOM voltages in Eq. I|33fl are therefore given 
by 

Vi(t)= ^Wiit-t'^WAt^-x^dt' . (67) 

J-oo I=1 

With the manipulations that led to Eq. I|37|l this leads to 

n(t,t[) N p 

Vi(t) = I Q t b Y,W l (t - t( - mt b ) (68) 

m— — oo 1=1 
I N L ,{i-l) 



X 



E E T^vMh+tj) 

J=l j=l 



x 



Of 



Following the derivation to Eq. (|44|l leads to 



N„ 



1=1 

I N L ,(i-l) 

e e n 

J=l j=l 



ii E 



iunt b 



n— — oo 



(.69) 



For all oscillation frequencies that are not multiples of the 
bunch repetition frequency, w = 2ivl/tb, I is an integer, 
the term xL Y]°°_ „ e lw ™*t vanishes so that the condition 
for Vt(oj) to be non-zero is the same as for the BBU 
instability without misalignments x^. Below threshold, 
the HOM voltage therefore has the following form: 



VM = £ v ie l 



(70) 



The voltage seen by bunch i on turn L is therefore 



given by V t L = £i°!Loc ^ ' 
t\ , Eq. JHEJ can be written as 



so that for the time 



vt = i tb 

771 — — OO 1—1 

I Njj(i-1) 

E E 

J=l 3 = 1 



t\ - mtb) (71) 



T IJ -V 



J 



' Qi 



After performing the summation over m one obtains 

N p 

Vt = / ^<=°(^,L)) (72) 
i=i 

I N u (i-1) 

E E >-!:!^' •<•• 

J=l 3 = 1 

where the superscript u> — means that w(8) is computed 
with Eq. (O for u = 0. 

Comparing with Eqs. (|44|> and 1)45(1 shows that this can 
be written as 



v = / w(o)uy - i w(o)x , 



(73) 



with W = W(0) and U from Eqs. gJJ and gSJ. The 

vector of voltages is therefore given by V = [JqWU — 
I] _1 /oW.To- Equation l-il'li can now be used to compute 
the beams distance from the cavity center at each turn 
and each cavity. One obtains 



/ N u (i-1) 



rplj Ty a 



b 0i •> 



(74) 



4ti)-4i = E E 

J=l 3 = 1 ° 

U[/ WU - IJ^/oWfo - x (75) 



X — Xq 



[JoUW-I]- 1 ^ 



Evaluating this for a single cavity with a single HOM 
leads to 



X — Xq 



•i'O 



f I t b T 12 w"=°(5) - 1 



(76) 



For S — h and with e = ^-tb one obtains 



X — Xq = Xq 



hKTv. 



cosh ; 



coshe — cos uj\tb 



(77) 



l=— oo 



There is a current Iq at which the denominator be- 
comes and the orbit deviation would become very large. 
The question arises whether this current is larger than 
the BBU threshold 7 t h or smaller, so that large orbit ex- 
cursions would present a new kind of instability. 

This problem does not only arise for the single HOM 
case of Eq. (|75|l but also for the general case of Eq. I|77|l . 
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Very large orbit excursions x occur for currents for which 
the matrix inverse does not exist. 

The inverse matrix to be inverted is (/jWU — = 
ioADiag[(Ai — l/Jo) -1 ] A -1 , where A is the matrix that 
diagonalizes W(0)U. We therefore see that for 
which the orbit gets very large is given by the eigenval- 
ues of W(0)U. These values are naturally smaller than 
1/ith, which is the largest eigenvalue of W(oj)U that is 
produced for any frequency to. 

This proves that the BBU instability always occurs 
before the orbit excursion becomes very large. 

VIII. TRACKING RESULTS 

The tracking code BI (stands for beam instability) was 
developed to perform studies of beam breakup in recir- 
culating linacs [2j|. The algorithm models point charge 
bunch interactions with HOMs in linacs, taking into ac- 
count proper time delays between the cavities, transfer 
maps, etc., allowing BBU simulations due to longitudi- 
nal, transverse and other higher order modes in a general 
linac configuration. 

The basic algorithm can be summarized as following. 
The string of HOMs that a bunch sees in its lifetime be- 
tween injection and ejection points is represented by a list 
of pointers to the actual cavities. The proper time delays 
between cavities is also stored for each pointer. E.g. for N 
HOMs and N p passes, the list of pointers would be NN p 
long pointing to N HOMs. This approach allows one to 
represent any recirculation configuration without limita- 
tions. As the train of bunches is injected into the struc- 
ture, the next instance when any bunch sees any pointer 
is determined, and the HOM voltage in the correspond- 
ing cavity is updated. Then, this bunch is pushed to the 
next pointer where its coordinates are stored, waiting 
for its turn in time to be the next bunch going through 
a pointer. This way no bunches end up ahead of time 
precluding a situation when a bunch sees a cavity with 
incorrectly updated HOM fields, i.e. causality is properly 
realized. Furthermore, the algorithm is general enough to 
allow modeling of the longitudinal instability where tim- 
ing between different bunches is no longer kept fixed. The 
practical realization of this algorithm is relatively fast, al- 
lowing the tracking of a complete 5 GeV ERL with 300 
HOMs for 0.1 ms in less than a minute on an average per- 



sonal computer. This duration is sufficient to determine 
the onset of transverse BBU instability in most practical 
cases. 

The output of the code contains amplitudes of HOM 
voltages as a function of time, which is used to determine 
the growth rate of the instability by fitting an exponen- 
tial. Several successive calls are made to the tracking 
unit to determine the threshold. The length of the track- 
ing time is estimated from Eq. I|30l) based on the desired 
accuracy in threshold determination. 

IX. CONCLUSION 

For dipole HOMs we have derived the BBU theory for 
arbitrary recirculation times, so that the theory can be 
applied to ERLs. The resulting equations have been used 
to find analytical results for (1) multiple HOMs in one 
cavity, (2) two equal HOMs for a one-turn ERL, and 
(3) one HOM for a multiple times recirculating ERL. 
For (1) the numerical observation that it often suffices to 
include only the strongest of several different HOM was 
explained and the distance in frequency was derived for 
which one can consider two HOMs as different, for (2) 
it was shown that two cavities do not cancel each others 
instability but that they can be arranged so that they 
do not add dangerously, and for (3) was shown that the 
BBU threshold current for an N r times recirculating ERL 
is roughly up to N r (2N r — 1) times smaller than that in 
a corresponding one-turn ERL. 

Furthermore, a simple method to compute the orbit de- 
viations produced by cavity misalignments has also been 
introduced. And it is shown that the BBU instability 
always occurs before the orbit excursion becomes very 
large. 

Several comparisons with tracking data verify the ap- 
plicability of the theory and of the tracking program. The 
conclusions should be useful in determining and optimiz- 
ing the maximum current for the many ERLs that are 
currently in design and pre-proposal stages worldwide. 
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